Draft version February 1, 2008 

Preprint typeset using liT^ style emulateapj v. 04/03/99 



O 
O 

o 

> 

o 
en 

> 

m 

o 
o 

Oh' 
I ' 



% 



IMPLICATIONS FOR QUINTESSENCE MODELS FROM Maxima-1 AND BooMERANG-98 

A. Balbi', C. BACciGALUPr, S. Matarrese'''*, F. Perrotta", N. Vittorio' 
Draft version February 1, 2008 

ABSTRACT 

Prompted by the recent Maxima-1 and BoOMERANG-98 measurements of the cosmic microwave background 
(CMB) anisotropy power spectrum, and motivated by the results from the observation of high-redshift Type la 
supernovae, we investigate CMB anisotropics in quintessence models in order to characterize the nature of the dark 
energy today. We perform a Bayesian likelihood analysis, using the Maxima-1 and BoOMERANG-98 published 
bandpowers, in combination with COBE/DMR, to explore the space of quintessence parameters: the quintessence 
energy density ft^ and equation of state w^. We restrict our analysis to flat, scale-invariant, inflationary adiabatic 
models. We find that this simple class of inflationary models, with a quintessence component 17^ < 0.7, -1 < 
W0 < -0.5, is in good agreement with the data. Within the assumptions of our analysis, pure quintessence models 
seem to be slightly favored, although the simple cosmological constant scenario is consistent with the data. 

Subject headings: Cosmic microwave background — Cosmology: observations — Cosmology: theory — 
equation of state 



1. INTRODUCTION 

Recent results (Perlmutter et al. 1999; Riess et al. 1998) 
from the observation of high-redshift Type la supernovae (SNe 
la) have driven a renewed interest towards cosmological sce- 
narios where the bulk of the energy density is made up of some 
unknown negative-pressure component (named "dark energy" 
in the recent literature). Among the possibilities that have been 
proposed to explain the nature of this component is the exis- 
tence of a non-zero cosmological constant (constant vacuum 
energy). A more general way to describe this component is 
by its equation of state p = wp, with w < 0. This allows for a 
whole set of candidates with different values of w. Among them 
is a dynamical form of time varying, spatially inhomogeneous 
dark energy component, generally provided by a slowly rolling 
scalar field, referred to as "quintessence" (Caldwell et al. 1998). 
The most attractive feature of quintessence models is their ca- 
pability to solve the so-called 'cosmic coincidence' problem, 
namely the fact that matter and vacuum energy contribute by 
comparable amounts to the energy density today. Successful 
quintessence models should in fact allow for the existence of at- 
tractor (Ratra & Peebles 1988) or 'tracker' (Zlatev et al. 1999) 
solutions, which make the present-day behavior of the scalar 
field nearly independent of its initial conditions (e.g. Wang et 
al. 2000, for a thorough review). 

It has been noticed (see, e.g., Perlmutter, Turner & White 
1999) that it may be difficult for SNe la data alone to discrim- 
inate among different candidates for the dark energy. Comple- 
menting these data with those coming from other observations 
may help assessing the nature of the dark energy component. 
Study of the cosmic microwave background (CMB) anisotropy 
has long been recognized as a powerful tool to investigate cos- 
mological models. The recent outstanding achievements of the 
Maxima-1 and Boomerang-98 experiments (Hanany et al. 
2000; de Bernardis et al. 2000) have finally moved us closer 
to the goal of a long awaited high precision measurement of the 
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CMB angular power spectrum. The CMB power spectrum mea- 
surements by Maxima-1 and Boomer ang-98 are in remark- 
able agreement, and have akeady been used to set constraints 
on the values of a suite of cosmological parameters within the 
class of inflationary adiabatic models (Jaffe et al. 2000; Balbi 
et al. 2000; Lange et al. 2000). These results, which strongly 
support a universe with a density very close to critical, made 
the case for the existence of dark energy even stronger. In 
fact, large scale structure observations indicate that matter with 
strong clustering properties can only make up to ^ 30% of the 
critical density, leaving room for a ^ 70% contribution from 
dark energy. It is then timely to use these new CMB datasets 
to explore in more detail the nature of the dark energy. A com- 
parison between quintessence models and recent CMB data has 
been performed by Brax et al. (2000); they investigated the de- 
pendence of CMB power spectrum (in particular the location 
and height of acoustic features) on the main cosmological pa- 
rameters as well as on different quintessence potentials. The 
predicted CMB anisotropy is then compared to BOOMERANG- 
98 and MAXIMA- 1, finding general agreement. A likelihood 
analysis of the BoOMERANG-98 results for theories involv- 
ing couplings of the quintessence field with other cosmological 
components has been carried on by Amendola (2000), finding 
again broad compatibility with the data. In this Letter we focus 
on the quintessence scenario, and use the CMB power spec- 
trum from Maxima-1 and Boomer ang-98 to set Hmits to 
the parameters of this model. These are the present-day clo- 
sure energy density of the quintessence component, fi^, and 
the quantity w^ = p^j p^, characterizing its equation of state. 
We have purposely restricted our analysis to the simplest class 
of inflationary theories, namely flat models with scale-invariant 
adiabatic scalar perturbations. 

This Letter is organized as follows. In Section 2 we give a 
quick review of the key features of the quintessence scenario, 
with a particular emphasis to the theoretical predictions for 
the CMB angular power spectrum. In Section 3 we use the 



Maxima- 1 and Boomerang-98 power spectrum measure- 
ments to constrain the parameters of the model. Finally, we 
discuss our results in Section 4. 

2. CMB ANISOTROPY IN QUINTESSENCE MODELS 

The quintessence component is usually described as an ultra- 
light (its Compton wavelength being of the same order or larger 
than the Hubble radius) self-interacting scalar field (p, with neg- 
ative pressure. It is invoked to provide the accelerated expan- 
sion of the Universe today in alternative to the ordinary cosmo- 
logical constant. In this scenario, the vacuum energy is stored 
in the potential energy V((j>) of the (p field, evolving through 
the classical Klein-Gordon equation. A suitable period of slow- 
rolling in the scalar field dynamics allows to reproduce a behav- 
ior similar to a cosmological constant, which in turn requires 
that the potential should admit a rather flat region. Several 
candidates have been proposed in recent years in the particle 
physics context (e.g. Masiero, Pietroni & Rosati 2000, and refs. 
therein). In this Letter we focus on quintessence models with 
inverse power-law potentials, V oc (/)"" (a > 0), as originally 
suggested by Ratra & Peebles (1988). 

The large interest in quintessence models is due to their ca- 
pability to avoid the severe fine-tuning on the smallness of the 
vacuum energy compared with matter and radiation in the early 
universe: a large class of field trajectories, named "tracking so- 
lutions" (Zlatev et al. 1999), are able to set the energy at 
the observed value today starting from a very wide set of ini- 
tial conditions; among others, Ratra-Peebles-type potentials do 
satisfy this requirement (Liddle & Scherrer 1999). 

For a given class of potentials, quintessence models can be 
parametrized by two quantities: the scalar field closure energy 
density today, 17^, and its equation of state at the present time, 
W0 = [p<f,/P4,]to^ where pi/, and p^ are the field pressure and en- 
ergy density, respectively. Notice that w^ = -l for a pure cos- 
mological constant, but can be larger for scalar field compo- 
nent. In particular, in the tracking regime p^ ex (1 +2)3"/("+2) ^^ 
low redshifts, thus reducing to the ordinary cosmological con- 
stant for a ^ 0. Moreover, there is a one-to-one relation be- 
tween Wcf, and the potential exponent a: w^ = -2/(a + 2): the 
larger w^, the greater the inverse potential exponent a. How- 
ever, one must be aware that, although useful for a qualitative 
understanding, the above relations are only approximately sat- 
isfied at present, since the quintessence field is dominating and 
the tracking regime is being abandoned. 

Even though the (j) field is too relativistic to cluster on sub- 
horizon scales, its dynamics induces relevant imprints on cos- 
mological scales. The scales beyond the horizon at decoupling 
are depicted on the \ow-l tail of the CMB anisotropy angu- 
lar power spectrum; the quintessence dynamics acts on these 
scales through the Integrated Sachs-Wolfe effect (ISW; e.g. Hu, 
Sugiyama & Silk 1997). Also, the (j) component participates in 
the overall expansion of the Universe, causing a geometric ef- 
fect that translates into a projection of each physical scale onto 
different angles in the sky. Moreover, since its energy density 
was largely subdominant at decoupling, the form of the acoustic 
peaks region of the CMB spectrum at £ > 200 is not altered. 

For a given value of the vacuum energy density today, the 
effects of a quintessence component, compared to the cosmo- 
logical constant, can be qualitatively understood by consider- 
ing the dark energy equation of state. When moving away from 
W0 = -1, the first effect that is important is the projection one: 
the size of the last scattering surface gets reduced because in 
the past the quintessence energy density was larger than that of 



a cosmological constant: this determines a shift of all the CMB 
features to larger scales, or smaller multipoles. In addition, the 
more w^ deviates from -1, the earlier is the time when (f) starts 
to dominate the cosmic expansion; through the ISW effect, the 
CMB photons are sensitive to the dynamics of the gravitational 
potential wells and hills along their trajectories, which in turn 
decay much more rapidly, reflecting the dominance of the vac- 
uum energy. This enhances the level of CMB anisotropics on 
the low £'s tail of the spectrum, where one should consider the 
limitations introduced by cosmic variance, when dealing with 
the ISW power However, for the tracking solutions analysed 
here, the variation in the power spectrum introduced by the ISW 
for -0.6 <w^< -0.5 is larger than the uncertainty due to cos- 
mic variance. 

CMB anisotropy angular power spectra in the considered 
class of quintessence models have been obtained by suitably 
modifying CMB FAST (Seljak & Zaldarriaga 1996) to include 
(f) perturbations in the synchronous gauge, as detailed in Per- 
rotta & Baccigalupi (1999), and numerically evolving the back- 
ground quantities, accounting for the tracker scalar field tra- 
jectories, as discussed in Baccigalupi, Matarrese & Perrotta 
(2000), where results were presented for some specific values 
of the quintessence parameters. 
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Fig. 1. — Data points used in the analysis. The full circles 
are the bandpowers from MAXIMA- 1, the open circles are from 
BoOMERANG-98. These data have been re-scaled within their \a cal- 
ibration uncertainties. The open boxes are the COBE/DMR data. The 
curves are the theoretical predictions from quintessence models which 
best fit the combined data when we allow the physical baryon density 
to vary in the range 0.01 < Q,tl'i^ < 0.05 (continuous line, correspond- 
ing to Q,^ = 0.5, W0 = -0.75 and Q.bh'' = 0.03) or we fix it at the big 
bang nucleosynthesis value Q.bh^ = 0.02 (dashed line, corresponding 
to Q,^ = 0.7 and w^ = -0.75). 

3. COMPARISON WITH THE DATA 

In this section we compare the predictions of quintessence 
models to the results of the first flight of the MAXIMA ex- 
periment (Maxima- 1; Hanany et al. 2000) and of the long 
duration Antarctic flight of BOOMERANG (BoOMERANG-98; 
de Bernardis et al. 2000). To rule out models in the 2- 
dimensional space of quintessence parameters (ft^, w^), we 
adopted a Bayesian approach, estimating the probability dis- 
tribution of the parameters by evaluating the likelihood of the 



data, L. oc exp(-x^/2). The data we used are the MAXIMA- 1 
spectrum, estimated in 10 bins covering the range 36 < i? < 785, 
and the BoOMERANG-98 spectrum, estimated in 12 bins, over 
the range 26 < i' < 625. Within each bin, the power spectrum 
is parameterized by a flat bandpower: C/, = i{l+ 1)Q. We used 
the published bandpowers and neglected the effect of bin-bin 
correlations (which should be small, see Hanany et al. 2000, 
de Bernardis et al. 2000), and the offset-lognormal corrections 
to the likelihood, since these have not been made public yet. 
Foflowing Jaffe et al. (2000), we combined the MAXIMA- 1 
and BoOMERANG-98 adjusting the calibration of each dataset 
by a factor 0.98 ±0.08 for Maxima-1 and 1.14 ±0.10 for 
BooMERANG-98. We extended the low €•& coverage of the 
dataset by including in our analysis the COBE/DMR power 
spectrum estimated over 8 bins in the range 2 < i? < 22 by 
Tegmark & Hamilton (1997). The band powers used in the 
analysis are shown in Figure 1 . 

Since the aim of this work is to focus specifically on the 
quintessence parameters space, we did not embark in a multi- 
dimensional fit such as those already performed in Lange et 
al. (2000), Balbi et al. (2000), Jaffe et al. (2000). We inten- 
tionally decided to restrict our analysis to the most basic set 
of inflationary adiabatic models, namely those with a flat ge- 
ometry (il = fi^ + il,,, = 1) and a scale-invariant spectrum of 
primordial fluctuations (n., = 1). Furthermore this class of infla- 
tionary models is precisely the one which seems to be favored 
by the current multi-parameter CMB analyses (see also, e.g., 
Kinney, Melchiorri & Riotto, 2000). We treated the overall 
normalization of the power spectrum as a free parameter, and 
for each model we adjusted it to the value which gave the best 
fit to our data. We fixed the value of the Hubble constant at 
i/o = 65 km s"' Mpc"', which is consistent with most of the 
observations (Freedman 1999). We assumed that reionization 
occurred late enough to have a negligible effect on the power 
spectrum (corresponding to an optical depth r = 0). Finally, we 
did not include the effect of massive neutrinos (which would be 
very small anyway) and of tensor modes (gravitational waves). 

We only made an exception to our strategy of using very tight 
priors in the case of the physical baryon density il^/!^ (where h 
is the Hubble constant in units of 100 km s"' Mpc"'), since 
the value VL\,h^ = 0.019 ± 0.002 deduced from big bang nucle- 
osynthesis (BBN) considerations using the measured deuterium 
abundances (Buries et al. 1998) seems to conflict with the value 
nbh^ = 0.032 ±0.005 (Jaffe et al. 2000) emerging from the al- 
ready mentioned analyses of the CMB datasets (see also Es- 
posito et al. 2000). Then, in order to test the importance of this 
discrepancy in our analysis, we allowed the physical baryon 
density to vary in the range 0.01 < flhh^ < 0.05, which en- 
compasses a more conservative range consistent with high and 
low deuterium abundances measurements. In this case our re- 
sults were obtained by marginalizing (integrating) the likeli- 
hood over il.hh^. 

Our results are shown in Figure 2 where we plot the 68% 
and 95% confidence levels in the fJ^-w^ plane. When we let 
flhh^ vary, we are able to find quintessence models which are 
an excellent fit to our data. We point out that most of the weight 
to the likelihood is provided by models with flbh^ ~ 0.03. In- 
deed, marginalizing over ilhli^ or fixing it at rihh^ = 0.03 give 
almost exactly the same result. The best fit is a model with 
n^ = 0.5, W0 = -0.75, flhh^ = 0.03, having x^/DOF=26/27 
(where DOF indicates the degrees of freedom, resulting from 
30 data points and 3 parameters). If we fix the physical baryon 



density to the BBN value f7/,/z^ = 0.02, the bounds we obtain 
in the il^-w^ plane are qualitatively similar to those of, e.g., 
Perlmutter, Turner & White (1999), Efstathiou (1999), Wang et 
al. (2000). However, we stress the fact that none of the models 
with ilhli^ = 0.02 provides a good fit to the data. In fact, the best 
fit, with n^ = 0.7 and w^ = -0.75, has x^/DOF=38/27. QuaUta- 
tively, we found that increasing H^/z^ has the effect of moving 
the best fit to lower values of H,^, in a roughly linear fashion, 
while W0 remains roughly constant at w^ ^ -0.75. The best fit 
models with i^hh^ = 0.03 and fi^/z^ = 0.02 are compared to the 
data in Figure 1 . 
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Fig. 2. — Likelihood contours at 68% and 95% confidence level 
in the Q^ — w^ plane. The shaded regions are obtained marginaliz- 
ing over Qhh in the range 0.01 < flth < 0.05. For comparison, the 
dashed contours show the bounds obtained imposing the big bang nu- 
cleosynthesis value Qbh^ = 0.02. 
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Fig. 3. — Marginalized likelihood, normalized to the peak value, 
for the individual quintessence parameters Q^ and w^. In both pan- 
els, the solid line is obtained marginalizing over flhh^ in the range 
0.01 < flth^ < 0.05, while the dashed line is obtained imposing the 
big bang nucleosynthesis value Q,th^ = 0.02. Note, however, that the 
absolute likelihood obtained for Q,th^ = 0.02 would be more than two 
orders of magnitude lower than that obtained marginalizing over Q,i,h^ 
(see text). 



In Figure 3 we show the projected likelihood functions for 
the individual quintessence parameters, obtained by marginal- 
ization of the 2-dimensional likelihood with respect to each pa- 
rameter. A robust indication, quite independent of the value 
of ilhh^, is that although pure cosmological constant models 
(w0 = -1) are compatible with the CMB data used here, pure 
quintessence models (w^ > -1) seem to be slightly favored. 
This appears to agree with the qualitative discussion of Brax, 
Martin & Riazuelo (2000). The peak location in models with 
W0 > -1 is shifted to the left with respect to models with 
Wcf, = -1, resulting in better agreement with the data. In fact, 
this mechanism allows one to find quintessence models with 
J7 = 1 which have the same peaks locations of non-quintessence 
models with a il slightly larger than unity. If we impose 
W0 = -1, we get a best fit with il^ = 0.45 and Ui,h^ = 0.03 
with xVDOF=28/27 (n^ = 0.65 if we impose f^^/i^ = 0.02, 
with x^/DOF=40/27). It is also interesting to note that these 
data seem to disfavor high values of H.^, unless we impose the 
BBN constraint on i^bh^- However, we should keep in mind 
that this result is influenced by the assumptions we made on the 
value of the other parameters. Relaxing the priors would result 
in higher upper bounds on fJ^, similar to those found in Lange 
et al. (2000) or Jaffe et al (2000). 

4. DISCUSSION AND CONCLUSIONS 

We used the recent CMB power spectrum measure- 
ments from Maxima-1 and BoOMERANG-98, combined with 
COBE/DMR, to investigate the quintessence models in the con- 
text of flat, scale-invariant, inflationary adiabatic models. The 
Maxima- 1 and Boomerang-98 data already provide inter- 
esting insights on the equation of state of the universe. The sim- 



plest class of inflationary models, with quintessence, together 
with economic assumptions on the other parameters, provide 
an excellent fit to these CMB data. In fact, quintessence models 
seem to be slightly favored over simple cosmological constant 
models. 

In this Letter we restricted our attention to the simplest class 
of quintessence potentials, having the Ratra-Peebles V oc </)"" 
form. Within this class of models one can already conclude 
that shallower potentials (a < 2) are favored by current CMB 
anisotropy data. Also, we did not analyze here the case 
when the scalar field providing the dark energy today is non- 
minimally coupled to the gravitational sector of the theory, as it 
happens in the 'extended quintessence' scenario (Perrotta, Bac- 
cigalupi & Matarrese 2000; Baccigalupi et al. 2000). We did 
not compare explicitly the constraints on the equation of state 
coming from high-redshift SNe la observations (Perlmutter, 
Turner & White 1999; Garnavich et al. 1998) with the results 
of our analysis, because of our choice of considering a range of 
values for ri^/z^ not limited to the standard BBN value. Clearly, 
high precision SNe la distance measurements such as those ex- 
pected to come from the proposed SNAP^ satellite will provide 
powerful constraints on quintessence models complementary to 
those of the MAP* and Planck^ CMB missions. Similar com- 
plementary information will come from measurements of the 
comoving density of galaxy clusters (Haiman, Mohr & Holder 
2000). Thus, one will be able to put severe constraints not only 
on the nature of the dark energy, but also on the structure of the 
theory of gravity. 
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